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Abstract 

The singular parabolic problem Ut = Au — (^^.^^^y on a bounded domain Q. of with 
Dirichlet boundary conditions, models the dynamic deflection of an elastic membrane in 
a simple electrostatic Micro-Electromechanical System (MEMS) device. In this paper, we 
analyze and estimate the quenching time of the elastic membrane in terms of the applied 
voltage — represented here by A. As a byproduct, we prove that for sufficiently large A, finite- 
time quenching must occur near the maximum point of the varying dielectric permittivity 
profile f{x). 
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1 Introduction 

Micro-Electromechanical Systems (MEMS) are often used to combine electronics with micro- 
size mechanical devices in the design of various types of microscopic machinery. An overview 
of the physical phenomena of the mathematical models associated with the rapidly developing 
field of MEMS technology is given in [12]. The key component of many modern MEMS is the 
simple idealized electrostatic device shown in Figure [H The upper part of this device consists 
of a thin and deformable elastic membrane that is held fixed along its boundary and which 
lies above a rigid grounded plate. This elastic membrane is modeled as a dielectric with a 
small but finite thickness. The upper surface of the membrane is coated with a negligibly thin 
metallic conducting film. When a voltage V is applied to the conducting film, the thin dielectric 
membrane deflects towards the bottom plate, and when V is increased beyond a certain critical 
value V* -known as pull-in voltage- the steady-state of the elastic membrane is lost, and proceeds 
to quenching, i.e. snap through, at a finite time creating the so-called pull-in instability. 



'Partially supported by the Natural Science and Engineering Research Council of Canada. 



1 



EfelrMtiiHEwtCnliig 



SLfpafidBmiy 



KaljiQiiSte 





y 



Figure 1: The simple electrostatic MEMS device. 



A mathematical model of the physical phenomena, leading to a partial differential equation 
for the dimensionless dynamic deflection of the membrane, was derived and analyzed in [7]. In 
the damping-dominated limit, and using a narrow-gap asymptotic analysis, the dimensionless 
dynamic deflection u = u{x, t) of the membrane on a bounded domain 0, in R2, is found to 
satisfy the following parabolic problem 



ut — Au 

u{x, t) 
u{x, 0) 



A/(x) 



(1 







for X £ Q , 

for X E . 
for X £ il.. 



iP)> 



The initial condition in {P)x assumes that the membrane is initially undeflected and the voltage 
is suddenly applied to the upper surface of the membrane at time t = 0. The parameter A > in 
{P)x characterizes the relative strength of the electrostatic and mechanical forces in the system, 
and is given in terms of the applied voltage V hy X = ^"^J; , where d is the undeflected gap 
size, L is the length scale of the membrane, Tg is the tension of the membrane, and Eq is the 
permittivity of free space in the gap between the membrane and the bottom plate. We shall 
use from now on the parameter A and A* to represent the applied voltage V and pull-in voltage 
V*, respectively. Referred to as the permittivity profile, f{x) in {P)\ is defined by the ratio 
f{x) = ^^jj where £2(2;) is the dielectric permittivity of the thin membrane. 
Consider first the steady-state solutions of {P)x 



-Aw 

w^x) 



A/(x) 



w 



x£dn 



with 0<t(;<lonOc M , and /(x) was assumed to satisfy 

/ G C"(J^) for some a G (0, 1], < / < 1 and 
/ > on a subset of J7 with positive measure. 



2 



One can then easily show (e.g., Theorem 1.1 in [4j) that there exists a finite pull-in voltage 
A* := X*{n,f) > such that: 

• If < A < A*, there exists at least one solution for (S')a- 

• If A > A*, there is no solution for {S)x. 

Upper and lower bounds on the pull-in voltage A* were also given in Theorem 1.1 of [4J. Fine 
properties of steady states -such as regularity, stability, uniqueness, multiplicity, energy esti- 
mates and comparison results- were shown in ^ and [53 to depend on the dimension of the 
ambient space and on the permittivity profile. 

For the dynamic problem {P)x, we first define the following notion. 
Definition 1.1. (1) A solution u{x,t) of {P)\ is said to be quenching at a -possibly infinite- 
time T = T(A, /, il), if the maximal value of u reaches 1 at time T. 

(2) A point xq £ 0, is said to be a quenching point for a solution u{x, t) of {P)x, if for some 
T £ (0, -|-oo], we have lim u{xQ,tn) = 1. 

In [5] we dealt with issues of global convergence as well as quenching in finite or infinite time 
of the solutions of {P)x- One of the main results was the following relationship between the 
voltage A and the nature of the dynamic solution u of (-P)a- 

Theorem A (Theorem 1.1 in [5j). Assuming that f satisfies (1.1) on a bounded domain 0,, 
then the fallowings hold: 

1- If ^ ^ then there exists a unique solution u{x,t) for {P)x which globally converges 
pointwise as t ^ +oo to its unique minimal steady-state. 

2. If X > A* and inf^j / > 0, then the unique solution u{x, t) of {P)x nnust be quenching at a 
finite time. 

A refined description of finite-time quenching behavior for u was given in [6] , where some quench- 
ing estimates, quenching rates, as well as some information on the properties of quenching set 
-such as compactness, location and shape, were obtained. 

The first purpose of this paper is to prove -in Theorem 2.1- that quenching in finite-time 
occurs as soon as A > A*, which means that Theorem A. 2. above holds without the restriction 
info / > 0. On the other hand, we continue our search for optimal estimates on quenching times 
at voltages A > A*, since the latter translate into useful information on the operation speed of 
MEMS devices. Indeed, we established in Theorem 1.3 of [5], that if inf^jgrj f{x) > 0, then the 
following upper estimate for the quenching time holds for any A > A*: 



Tx{^J) < 



8(A + A*)2 r /A + 3A*\i/2 



3inf,ef,/(x)(A-A*)2(A + 3A*) 



V2A + 2A*/ 



:i.2) 



In this paper, we shall improve this estimate -at least in dimensions less than 8- by proving 
that 

Tx{nj)r^C{X-X*)--^ as A\A*, 
while ^ 

^ ~ T{ TTT A / oo. 

3Asup^gf^/(x) 
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To be more precise, we first recall the following notions and results from [4J . For any solution 
w of {S)x, we consider the linearized operator at w defined by L^j^x = — ^ — (i-w)^ ' ^^^^ 
corresponding eigenvalues {fJ'k,x{w);k = 1,2,...}. Say that a solution wx of {S)x is minimal, if 
wxix) < w{x) in whenever w is any solution of {S)x- We recall the following 
Theorem B (Theorem 1.2 in ||4j). Assume f satisfies (ji.ip on a bounded domain U C M^. 
Then, 

1. For any < A < A*, there exists a unique minimal solution wx of {S)\ such that Hi^x{w\) > 0. 
Moreover for each x G il., the function A — > wx{x) is strictly increasing and differentiable on 
(0,A*). 

2. If 1 < N < 7, then w* = lim wx exists in C^'^(Q) which is then a solution for (S)x* such 

\]\* 

that /Ui,A*(^^*) = 0. In particular, w* -often referred to as the extremal solution of problem {S)x- 
is unique. 

3. On the other hand, if N > 8, f{x) = \x\" with < a < a**{N) := 4-6jV+3y6(jV-2) ^ 
is the unit ball, then the extremal solution is necessarily w*{x) = 1 — |a;|~3~ and is therefore 
singular. 

We remark that in general, the function w* exists in any dimension, does solve {S)x* in a suitable 
weak sense and is the unique solution in an appropriate class. The above theorem says that it 
is however a classical solution in dimensions 1 < < 7, that is 



- Aw* = in 17, u;*>0 in 17 , w* = on dQ , (1.3) 



(1 - w*) 

and there exists an eigenfunction </>* of L^* x* satisfying 



2\*(h* f(x) 

A(p* + , \l =0 in 17 , (/>* > in 17 , = on 917 . (1.4) 

(1 — w*'Y 

We denote by (\)* (resp., V*) the corresponding unique L^-normalized (resp., L^-nornialized) 
positive eigenfunction of L^* x*. 

We shall then prove in section 2 the following upper and lower estimates on the quenching 
time T = r(A, /, 17) of a solution u for {P)x at voltage A > A*: Under the condition that the 
unique extremal solution w* of {S)x is regular, then 

• For A sufficiently close to A*, we have the lower bound estimate 



I. .-1 



T(A, /, 17) > (-— — ^ ""-gf^ l yi>^->^r-^- (1-5) 

If ^y^dx < oo, then for any A > A*, we have the upper bound estimate 

^ r 1 
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Note that the above situation typically happens when / = \x\^ and < 7, or for any N > 8 
provided /? is large. It would be interesting to establish similar estimates in the case where w* 
is singular. In the general case, we only have the following estimate established in section 3. 

• There exist a constant C = C{f,Q,) > and a sufficiently large Aq = Xo{f,Cl) > A* such 
that for any A > Aq, we have the estimates 

3XsnJ -f(x) - ^' ^) - 3Asup^- ffx) + ' ^^"^^ 

CSA SUp^gj^ / [X) SUp^gQ / {X) )^ 2+c 

where a E (0, 1] is as in (jl.ip . 

As a byproduct of the estimate p.7p . we shall analyze and compute in section 3 that in 
several situations, and at least for sufficiently large A, quenching in finite-time must occur 
near the maximum point of the varying dielectric permittivity profile /. More precisely, if the 
quenching set of a solution u for {P)\ is compact in Q, and if we are in one of the following 
two situations: 

1) N =1; or 

2) > 2, is a ball Bji{0), K = {0} and /(r) is radially symmetric, 

then for any a £ K, there exists C > such that for A large enough, we have 

(sup/)^-(/(a))^ <-^, (1.8) 

xf^Q A2+a 

We note that the compactness of the quenching set has been established in [6j (Proposition 2.1) 
in the case where the domain Q is convex and / satisfies both (jl.ip and the additional condition 

|i < on := {x G : dist{x, dn) < 5} for some 5 > 0. (1.9) 

Here u is the outward unit norm vector to dQ. The above result can be seen as a refinement of 
Theorem 1.1 of [6j where it is proved that under the compactness assumption on the quenching 
set, the latter set cannot contain any zero of the profile / (see also Lemma 3.2 below). 



2 Quenching time for A > A* 

In this section, we establish the estimates on the quenching time of (P)a- First we borrow 
ideas from [1] to prove that we have quenching in finite time as soon as A > A*, without the 
assumption used in [5] that / is bounded away from zero. 

Theorem 2.1. If X > X*{0,,f), then the unique solution u{x,t) of {P)\ must quench in finite 
time. 

Proof. The uniqueness of solutions for {P)\ in Q x (0, r), where r > is the maximal existence 
time, was already noted in Proposition 2.1 of [5j. Let now A > A*, and assume that u = u{x,t) 
of {P)x exists in X (0, oo). 
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Given any < e < A — A*, we first claim that {P)\-^ has a global solution that is uniformly 
bounded in x (0, oo) by some constant < 1. Indeed, set 

1 ds 
9iu) = ^^, hiu)=l ^, 0<n<l, (2.1) 

X — e ~ ds 
= In \i ' = / ^ ' < n < 1 , (2.2) 

and let <^e(u) := h^^ [h{u)) . Direct calculations show that 

-e A- 
■A + - 



$^(n) = 1 - + ^— ^(1 - n)^] ^ < Ce < 1 for 0<n<l 



Xfix)gi<^siu)) - {ve)t = \, ['y - ivs)t , 



1 

where Cg = 1 — (f) ^. Moreover, it is easy to check that ^^(O) = 0, that < ^e{s) < s for 
s > 0, and that $2(5) is increasing and concave with 

*'(.) = M>0. 

9{s) 

Setting Vs = ^e{u), we have 

-Ave = -#"(n)|Vn|2 - $',(n)Au 

> $^(^)(^M^ _ ^,) = Xf(x)<^'^{u)giu) - {Ve)t 

{X-e)f{_x) 

and hence, Vs = ^e{u) < Cg is therefore a supersolution of (P)A-e- Since now zero is a subso- 
lution of {P)x-£, we deduce that there exists a unique global solution for {P)x-£ satisfying 
< < < Ce < 1 uniformly in O x (0, 00), which gives our first claim. 
Note that (P)A-e admits a Liapunov functional 

V(.ue) = \ f \Vu,\^dx - (A - e) / -P^dx, V{ue) = - [ {u^ftdx. (2.3) 

Since now is uniformly bounded in $7 x (0, 00), we obtain that for /3 < 1, 

||(70,/3, ||uft||(^o,/3 < C uniformly bounded in X (0, 00) . (2-4) 

Moreover, (12. 3p gives that J^{ue)1dx < 00, which means that j^{ue)1dx is a uniformly 
continuous function on [0, co), and therefore 

{us)tdx as i — > 00 . 

Further, we deduce from ()2.4p that (ue)^ — > as i ^ 00, which shows that there exists a function 
< Ws{x) < Cg < 1 on Q such that Ue{x,t) W£{x) as t — > 00, where Ws satisfies 

-Au;, = ^^4?^ inO, we = on dQ . 

(1 - We)^ 

Therefore, there exists a classical solution of (5')a--£ with A — e > A*, which contradicts the 
definition of A* , and completes the proof of Theorem 12. 1[ ■ 
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2.1 Analytic estimates of quenching time 

We now focus on estimating the quenching time T when A > A*, and in the case where the 
unique extremal sokition w* of {S)\ is regular. This implies that w* satisfies 

- Au>* = / ^^^\, in 17 , w* >^ in 17 , w* = Q on dQ , (2.5) 

and there exists an eigenfunction (p* satisfying 
2X*(h* f(r) 

A^* + , \., =0 inQ, > inn, ^* = on dQ . (2.6) 

(1 — w*y^ 

We shall adapt and improve some of the arguments in [TO]. Our first estimate is a lower bound 
for T as stated in (1.5). 

Theorem 2.2. Suppose that the unique extremal solution w* of {S)\ is regular. Then for A 
sufficiently close to X* , the finite quenching time T(A, /, O) of the unique solution u for {P)x 
satisfies 

T{XJ,n)>( ^^P^^nr Y^x-X*)-K (2.7) 

H2A* sup.gf^ j^i^ ^ JT^dx^ 

where cp* > is the L'^ (Q) -normalized eigenfunction satisfying (2, 



Proof. Let u* be the unique solution of {P)x*- First, we seek a bound on the rate at which u* 
approaches the corresponding steady-state w* . For that, we set u*{x, t) = w*{x) — n(x, t). Then 
u{x,0) = w*{x) in 0, and u = w* on dQ. Moreover, we have 

du , . , , X*f(x) 
— = An - Aw* ■' ^ ' 



dt {l-w* + uf 

= Au + X*f{x)^ ^ ^ 

> Au + 



(1 — w*y {1 — w* + u)^ 

2X*uf{x) _ 3X*u'^f{x) 
(1 - w*)^ ~ (1 - w*)4 



(2. 



, . 2A*n/(2;) . 

> An + - Kiu^ , 

(1 — w*)-^ 



where Ki = 3A* sup^gj^ Define 



*=^, = (2.9) 

t + to -f^l 



where tg is chosen in such a way that 



V'(x,0) = — < n;*(j;) = n(x,0) in 0. 

^0 
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Note that I^M) gives 

and hence < ip < u = w* — u* in ^1 x (0, oo). 
We now set u = u* + ui, then ui satisfies 

du, , (A-A*)/(x) , , 
= Aui H h A f{x) 



dt (l-n)2 ^[(1-^)2 {1-u*) 

(A-A*)/(x) 2A*^i/(x) 
(l-u;*)2 ^ (l-u;*)3 ' 

as long as u = u* + ui < w* . We also define 

dx , F{x) = — < f(^x) 



(2.10) 



{l-w*f ' max{l,sup^gn/(x)} 
and consider <I>*(x) > to be a nonnegative solution of the problem 

+ (^*- + (1^ - = ° " (2.11) 

^>*(2;) = o xeao. 

Consider also the function 

^1 = (A- A*)(/i(/)*t + $*) in 17x(0,r), (2.12) 

where r > is arbitrary. Then ■0i(x,O) = (A — A*)$* > = ui{x,0) in 0, and ^i(x,t) = = 
wi(x,0) on 30. Moreover, since F{x) < 1 in 0, we obtain from (j2.10p and (j2.1ip that 

(V'l - ui)t - A{ipi - ui) 
= (A - A*)/i(/>* - (A - X*)IitA^* - (A - A*)A$* - {ui)t + Am 

> (A - A*)/i0*(x) - (A - X*)hr{x)F{x) + ^^IM^ii;, - ui) 

(1 — w*)-^ 

2A*/(x) 

in J7 X (0,r), as long as u = u* + ui < w* . Therefore, the maximum principle implies that 
ipi > ui as long as u = u* + ui < w* . 
We now obtain that 

u = u* + ui<w* -i^ + il;i = w* — + (A - A*)(/i0*t + $*). (2.13) 

t + to 

But the right-hand side of (j2.13p is no larger than w* , provided that 

> (A- A*)(/i</)*t + ^>*) in n, 

t + to 
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which is equivalent to 



K2>{X-X*){t + to){ht + A), where A = sup 



It requires 
which is 

where 



(A - X*)Iit' + (A - A*)(/ito + A)t -K2 + A{\ - \*)to < 0, 

-(A-A*)(/ito + ^) + VA 
- 2/i(A-A*) ' ^ ^ 



A := (A - X*f{Iito + Af + 4/1 (A - X*){K2 - Ato{X - A*)). 
For A sufficiently close to A*, (j2.14p can be satisfied if 



Note that Tl is given by 



1 



n=( ^-P^en^^x) )^(A-A*)-. 

Therefore, we conclude from ()2.13p that u < w* in 0, x (0,Tl]. This implies that the finite 
quenching time T oi u satisfies T >Tl, and the proof is complete. ■ 

We now establish the upper bound on T as stated in (1.6). 
Theorem 2.3. Suppose that the unique extremal solution w* of {S)x is regular, and that 
In ^f(x) '^^ ^ where ip* > is the L^{0,) -normalized eigenfunction satisfying h2. 6\) . Then for 
any A > A*, the finite quenching time T = T(A, /, il) of the unique solution u for {P)\ satisfies 

V Svr / Jn f(x) "-^ 
4 \X* f^'il)*{x)f{x)dx 

Proof. Setting u = w* + v, then we have 

dv . , (A-A*)/(x) , X*f{x) 
- Aw* + H -7: h 



T(A, /, n) < ^( ,. /"//ir,, ) ' (A - A*) . (2.15) 



at (i-u)2 [i-(uj* + v)]'^ 

2X*vf{x) (A-A*)/(x) 
= ^-+(1^+ (1-n)^ ^'-''^ 



.[l-{w*+v)Y {l-w*Y {l-w*Y 
Multiplying (j2.16p by tjj* and integrating over J7, we obtain 

1 1 2v 



A* / rfix) 



[1-{W*+V)]'^ (1-W*)2 {l-W*f 



dx, 
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where ()2.6p is applied. We next define 

E{t) = [ il;*vdx, E{0) = - [ ^P*w*dx = -Eq € (-1,0); 
Jn Jn 



h 



[ ^*{x)f{x)dx < 



r{x)f{x) 

(l-n)2 



dx , lo 



3A* 



Using the inequahties 



1 



2v 

[1 - {w* +v)Y " (1 - w*Y ~ (1 - w*Y 
the Holder inequality yields that 



> 



7T if w > 0; 



1 



2v 



dx 



> 3X* / v^'il^*{x)f{x)dx > 



3X* 



n 

It follows from the above that 
dE 



Jn fix) "-^ 



'il;*vdx) =l2E^{t). 



dt 



> {X - X*)Ii + hE', S(0) = -Eo G (-1,0). 



We now compare E{t) with the solution F{t) of 



dF 



^ = iX- X*)Ii + hF^ , F(0) = -Eo G (-1, 0). 



(2.17) 



(2.18) 



Standard comparison principle yields that E{t) > F{t) on their domains of existence. Therefore, 

sup V > E{t) > F{t) . (2.19) 
n 

It is easy to see from ()2.18p that the quenching time Ti for F{t) is given by 
fi = (^ + arctan , ^l... ) ( (A - A*)/i/2' 



(A-A*)/i 



< 



n fix) 



dx 



4 \X* f^'ip*{x)f{x)dx 



(A -A*)- 



Therefore, for any A > A* the unique solution u of {P)x must quench at a finite time T 
T(A, /, O) < Ti, and we are done. 
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3 Quenching behavior for sufficiently large A 

In this section we discuss the quenching behavior of solutions of {P)x for A large enough. We 
begin with the following refined estimates for the quenching time as stated in (1.7). 

Lemma 3.1. Assume f satisfies (li. on a bounded domain Q, and suppose u is a quenching 
solution of [P)x at finite time T . Then, there exist a constant C = C{f,Q) > and a sufficiently 
large Aq = \o{f,0,) > such that for any A > Aq, we have 

<T< ^^_ _^ + __, (3.1) 



3Asup3.gn/(a;) SAsup^gf^/ 
where a G (0, 1] is as in (jj.jp . 

Proof. In order to obtain the lower bound of finite time T, we consider the initial value problem: 

dr]{t) _ AM 

~ir - {l-v{tW ' (3.2) 
77(0) = 0, 

where M = sup^-gQ /(x). From (j3.2p one has fo^^\^ ~ s)'^ds = t . If is the time where 
lim^^T, vi^) = 1) then we have T* = Jq{1 — s)^ds = Obviously, r]{t) is now a super- 

solution of u{x, t) near quenching, and thus we have 

1 1 



3AM 3Asup^.gf2/(x) ' 



which is true for any A > 0. 

We next prove the upper bound in (j3.ip . Let a G be such that /(a) = sup^g^ f{x), and 
suppose K = K{f, Q) is the Holder constant of /. Since / G C'^{Q) for some a G (0, 1], then for 
any sufficiently small e > 0, there exists 6 = (21?)^'^" such that 

f{x)>f{a)-^-, yxeQ:=B{a,6)nn, 

where B{d, 6) is a ball centered at a with radius 6. Let v be the solution of 

^t-^v= ^^ll''}'^ inQx(0,r.), (3 3) 

v{x,0)=0 in Q, v{x,t)=0 on dQ x {0,Ty ) , 



where Ty is the maximal existence time of (|3.3p . Comparison arguments shows that u > t; in 
Q X (0, Tm), where = min{T, T^}. Therefore, we have T <Ty. 

Our goal now is to estimate for sufficiently large values of A. Let ^i(5) be the first 
eigenvalue of —A in B(a,6), and let (p be the corresponding positive eigenfunction normalized 
such that Jq (j)dx = 1. Multiplying (j3.3p by and integrating over Q, we obtain 

[ <f>vdx= [ 4>Avdx + X{f{a) - ^) [ ^ dx 
dtjq Jq 2'jQ{l-vy 



-Hi{S) I (pvdx + X{f{a) - ^) I ^ x2 



(3.4) 

- d/r. 

{l-vf 
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Next, we define an energy-like quantity by E{t) = Jq (p^vdx so that £'(0) = and 

E{t) = I (j)^v dx < snpv / (pdx = su'pv. (3-5) 
Jo Q Jo Q 



Q 

Then, using Jensen's inequahty on the right-hand side of (|3.4p . we obtain 



Recall that there exists a constant D = D{N) > 0, depending only on N, such that fj,i{5) = 
D6~^. We now choose e = e(A, /, > such that 

a 2 

F 2 \ '2Jl2+a K 2+a 

^,(5)=Dd-' = D{—)-^ = ^e, i.e., e = . (3.6) 

Then there exists a sufficiently large Aq = Xo{f,f!,) > X* such that for any A > Aq, we have 
/(a) — e > and 

dE ^ A(/(a)-£) ^ Xe 

-dt- il-EY ^W^''^'^^ 



X{f{d)-e) Xe X{f{a)-e) 

This implies a finite quenching time Tg of E satisfying 

^1 1 C 



3A(/(a)-£) - 3A/(a) ' 
where C = C{f, is independent of A in view of ()3.6p . Therefore, we conclude from (13. 5p that 

1 C 

T <Ty <Te < - . H 2+2^ > 

3A/(a) _>^^+^ 

and the lemma is proved. ■ 
We now recall the following result proved in Theorem 1.1 of [B]. 

Lemma 3.2. Assume f satisfies (ji.jp for some a G (0, 1] on a bounded domain Q, C M^, and let 
u be a quenching solution of {P)x at finite time T . Assuming the quenching set of u is compact 
in Q, then 

1. No point a £ Cl satisfying f{a) = can be a quenching point of u; 

2. There exists a constant Af > such that 

M{T -t)^ <l-u{x,t) in nx{0,T). (3.7) 

The following result can now be seen as a converse of Lemma 3.2: for sufficiently large A, 
finite-time quenching must occur near the maximum point of the varying dielectric permittivity 
profile /. 
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Theorem 3.3. Assume f satisfies (ji.jp for some a G (0, 1] on a bounded domain 17 C M^, and 
suppose that u is a quenching solution of (-P)a <it finite time T , in such a way that the quenching 
set K of u is compact in fi. Then, for any a ^ K, there exists C > such that for A large 
enough, we have 

(sup/)^-(/(a))^ (3.8) 

provided we are in one the following two situations: 

1) N = 1; or 

2) N >2 and a = Vt is a hall B^{0) and f{r) is radially symmetric. 

Proof. The idea of the proof -inspired by [2]~ is to combine the estimates on quenching time 
given by Lemma l3.ll with the local energy estimates near any quenching point established in 
[6]. Given a quenching point a of n and its corresponding quenching time T, we define 

y = ^jT^ ^ s = -log(l-^), l-u{x,t) = {T -t)^w{y,s), 

then w satisfies 



pWs = V ■ (pVw) + -pw — 1 in n{s) X (0, oo) , 



1 Xpf{a + y T2e-r 

-pw 2 

3 w"^ 

where p{y) = e"'^'^/^ and Q{s) = {y : a + yT^e~^ G O}. The compactness assumption on the 
quenching set implies that there exists a sufficiently large sq > such that Bs{a) C il.{s) for any 
s > So- 

Consider now the "frozen" energy functional 

E{w) = U p\Vw\'dy-U pwHy-f M^dy , 

which is defined in the compact set Bs of ria(s) for s > sq. Note from Lemma [3 . 2 1 that /(a) > 0. 
Using the same argument of Lemma 2.10 in [6j, one can obtain 

[ p\ws\'^dy < [ pWs^dS + ^ [ p\Vw\'^{y ■ iy)dS 

Jbs ds Jqb^ dv 2s JdBs 



, 1 



Xpwsifia)- f{a + yT2e-2)] 

dE 

where 



dy (3.9) 

+ /1 + /2 + /3, 

as 



To estimate I3, we use Lemma 13.21 to infer that w has a lower bound, and since / € C"(r2), we 
apply Holder's inequality to deduce that 

h < CT^e-^' I p\yrwsdy < C7rfe-t^( / p\ws\^dyy . 

J Bs ^ Bs 
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Therefore, (|3.9p gives for s S> 1, 



dE 



ds 



.2 



^<-/ p\ws\''dy + CT^e~^'{ I p\ws\'^dy)' + Cs^e-'rH . (3.10) 



Maximizing now the right hand side of (|3.1U|) with respect to / p\ws\'^dy, it yields that for 

Jbs 

s > 1 

ds 

This leads to 

E{w) < E{w{y,0)) + 911 = E{T-\) + 



a a 

Under the compactness assumption on the quenching set, a proof similar to Theorem 1.3 in [6] 
(see also O [9] ) gives that 

lim w{y,s) = (3A/(a))^ := k{a) 

s— >oo ' 

uniformly on \y\ < C for any bounded constant C, and E{w{-, s)) E{k{a)) as s ^ cxd, provided 
one of the following conditions holds: 

1) N = 1; or 

2) N > 2 and a = 0, = Br{0) is a bounded ball and /(r) = /(|x|) is radially symmetric. 
Therefore, under the assumption of Theorem 13.31 we have the following upper bound 

E{k{a)) <E{T-z) + - . (3.11) 

a 

Observe that if 6 is a constant then the energy E can be rewritten as E{h) = TF{b), where 
r = J p{y)dy and F is the function 

F{z) = --z^ - , z>0. 
b z 

Since F attains a unique maximum at k{a) and F"(k{a)) = — 1, there exist 7 and P such that if 
\z - k{a)\ < 7 then F"{z) < -i, and if \F{z) - F{k{a))\ < f3 then \z - k{a)\ < 7. So we obtain 
from ([3lT]) that 

Fik{a)) < F(r"3) + 

a 

Choose Ai such that = {3. Then for A > max{Ao,Ai}, where Aq is as in Lemma l3.ll we 
have 

/3>- > F{k{a)) - F{T~3) . 

a 

Hence from the properties of F, we have k(a) — T~3 < 7, which implies F"{k{a)) < — ^. It now 
deduces from p. lip that 

> — > F{kia)) - F(T-^) > -[r-3 - k{a)f, 



aA" ~ a ~ ' ~ 4' 
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where Lemma l3.1l is applied in the first inequahty. This further gives that 

i C 



r-^-(3A/(a))^ <-^. (3.12) 

A2 



On the other hand, since Lemma |3. II gives 



^ 1 C I f C 

T<- r + ^X^ < TTT . 1 + 



3Asup^gf^/ y2+§- 3Asup^gj^/V \2+ 
we have 

C 



T-l > (3Asup/(x))3(l 



Therefore, we finally conclude that 

(sup/(x))^ - (/(a))^ < + ^ < ^ 

xen A3 + 2 A2+" A2+Q 

This completes the proof of Theorem [ 



Figure 2: Left figure (a): plots of 1 — u versus x at different times, where A = 10. Right figure 
(b): plots of 1 — u versus x at different times, where A = 100. 



Before ending this section, we now present a few numerical simulations on Lemma 3.1 and 
Theorem 13.31 Here we apply the implicit Crank-Nicholson scheme (see §3.2 of [7] for details), 
with the meshpoints = 6000, to (-P)a in the symmetric slab domain —1/2 < x < 1/2. We 
choose the varying dielectric permittivity profile f{x) satisfying 

r 1 - 16(x + 1/4)2 ^ if X < -1/4; 
f[a]{x) = l |sin(27rx)|, if |x|<l/4; (3.13) 

( 1 - 16(x - 1/4)2 ^ X > 1/4. 

Note that x = ±0.25 are two maximum points of f{x), and all assumptions of Lemma 3.1 and 
Theorem 13.31 are satisfied in view of (1.9). 
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Figure 3: Left figure (a); plots of 1 — u versus x at different times. Right figure (b): local 
amplified plots of (a) . 



Simulation 1. Quenching behavior for small A > A*: 

In Fig. 2(a): 1 — u versus x is plotted at different times for (-P)a at A = 10, where the quenching 
time is T = 0.05174132. The quenching is observed at x = ±0.204, a bit far away from the 
maximum points of profile f{x). In Fig. 2(b): 1 — u versus x is plotted at different times for 
{P)\ at A = 100, where the quenching time is T = 0.003523908. In this case, the quenching is 
observed at x = ±0.2535, very close to the maximum points of profile f{x). This simulation 
shows the necessary of the assumption that Lemma 3.1 and Theorem l3.3l hold only for sufficiently 
large A. 

Simulation 2: Quenching behavior for sufficiently large A: 

In Fig. 3(a), 1 — u versus x is plotted at different times for (P)\ at A = 10^, where the quenching 
time is T = 0.000003332783. In this case, two quenching points are observed at x = ±0.250165, 
more close to the maximum points of profile f{x). In Fig. 3(6) we show the local amplified 
plots of (a) near the maximum point x = 0.25 of f{x). By further increasing the value of A, we 
observe that quenching points become further close to the maximum points of f{x). 
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